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SUMMARY a 


Storage and signal limited mathematical models for the insulin 
secretion process are described. The use of mathematical models in 


research is discussed. Kinetic in vitro perfused pancreas studies using 


glucose as a stimulator of ‘insulin secretion are evaluated in terms of 


\ 


the models: The kinetics of secretion inhibition by diazoxide and 
diphenylhydantoin (DPH) are explained in terms of the models and the 
compartment model is modified to allow simulation of an "off response”. 
Essential constraints for a satisfactory model are described, The 
above results are discussed as they relate to observed beta cell 


electrical activity, ion fluxes and metabolite requirements. The dis- 


tribution and metabolism of insulin in vivo is studied experimentally 


and evaluated with models. These studies lead towards eventual syn- 


thesis of secretion models with ion flux models and whole body insulin 


kinetics models. 
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Abstract 


Storage and signal limited mathematical models for the insulin 
secretion process are described. The use of mathematical models in 
research is discussed. Kinetic in vitro perfused pancreas studies using | | 
glucose as a stinisint of insulin secretion are evaluated in terms of 
the models. The kinetics of gecretion inhibition by diazoxide and 2 


diphenylhydantoin (DPH) are explained in terms of the models and the 


compartment model is modified to allow simulation of an "off response". 
Essential constraints for a satisfactory model are described. The 

7 pone results are discussed as they relate to observed beta cell 
electrical activity, ion fluxes and metabolite requirements. The 
distribution and ssiebieh 48 insulin in vivo is studied experimentally 


and evaluated with models. These studies lead towards eventual syn- 


thesis of secretion models with ion flux models and whole body insulin 

kinetics models. — conclusions are made: 

1. Present insulin must satisfy a complex series 

of events including fast, transient increased or decreased release, 

later phase, and hypersensitivity. | 

2. Specific signal limited secretion models work as well as specific 
storage limited models. 


3. In addition to previously determined constraints (10) an addition- 


al hypersensitivity factor 0(G,t) (or equivalent elaboration) is 


required for sufficient insulin secretion models. 


(a) Q(G,t) has faster kinetics than the provisionary factor, P(G,t). 
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(b) Release of catecholamines from nerve endings in the in vitro 
perfused pancreas does not account for the fast decay of 
hypersensitivity factor (Q(G,t))at | 
4. The in vitro pancreas is sensitive to the rate Of change of glucose 
concentration. | 
Je The aveion of diazoxide, including overshoot behavior when the drug is 
discontinued, can be simulated by a 3 compartment model; whereas 
a two-compartment model is sufficient to describe the action of | 
DPH. 
6. Linear compartment models of in vivo insulin ‘distribution and 
clearance are not sufficient to simulate both high concentration 
bolus and steady infusion experiments with a single set of viene 


for the constants 


7. Models with a saturable pathway for the degradation of insulin 
are sufficient to simulate the clearance from the blood of insulin 
presented by high concentration bolus and by steady infusion using 
the single set of values for the constants. 
8. No present insulin distribution and clearance model with a saturable 
pathway is uniquely defined by the present data. 


9. Measurements of insulin levels following both a high concentration 


insulin bolus and a short steady infusion are necessary to adequately 
define the insulin metabolic characteristics of an individual a | 


10. Ac yination of a storage limited model and a signal limited insulin 


secretion model could simulate all of the experiments presented in 


\ 


this thesis, but such a model would not be unique. 


UMI 
A 


| 
| 
| 
| 
| 
| 
| 
| 
| 

| 

KS 


UMI 


phenomena have been implicated as key processes in the mechanism of 


process of glucose stimulated insulin secretion. However, one result 


even though the data is not sufficient to determine the problem's unique 


| hormones are secreted from the islets of Langerhans found in the pan- 


Le Introduction 


The following is an analysis of the process of insulin secretion 


using secretion kinetics and other relevant data. Many different 


secretion. These phenomena are listed and an attempt is made to con- 


struct realistic models which could account for secretion in terms of 


such processes. The models are characterized as signal limited or 


storage limited. In so-doing it is hoped that these models will proyide 


a structure upon which further research can be based. ~— 
& 


This thesis does not provide the complete and final model for the 
of this study is a list of constraints for all sufficient models. The" é 


work here shows how models are used to define the limits of a problem, 


and satisfactory solution. 


Though the pancreas has an endocrine and an exocrine function, in 


this thesis, only the endocrine function is discussed. Several endocrine 


creas. They include somatostatin, glucagon and insulin. This thesis 


will directly discuss only the secretion of insulin although many of 


the conclusions relate to the secretion of other hormones. Approximately 
1% of pancreatic tissue is in the form of islets. 80 or 90% of the cells 
in the islets are insulin secreting beta cells. Insulin is the only hor- | 
mone which can directly lower blood sugar concentration; a defect in the 


secretion or action of insulin can lead to the disease diabetes mellitus. 
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Insulin is a polypeptide hormone, which after synthesis is packaged 
- the beta cell vee membrane bound vessicles somewhere in the Golgi 
region. Immediately after packaging the eianianiia of the vessicles are 
not opaque but later as the contents crystalize dark granules appear 
under the ish eis, The bulk of the insulin in the beta cell 


seems to be stored in this membrane bounded granular form. Secretion 


is thought to occur via a process called mevinnia. In exocytosis 

the mesibrane of a vesicle fuses with the cell ewibenne and then the mem- 
brane breaks at the fusion point and ‘the granule is released into the 
extracellular space. 

In vitro pancreas perfusion studies of insulin secretion rate asa. 
function of time during constant square wave infusion of glucose shows 

a multiphasic response (Fig. 3) with a fast early phase followed by a 
later slowly rising rate of taeehae secretion. The late phase rises to 
a steady state in approximately 60 minutes (1-4). Similar multiphasic 
insulin secretion’patterns during constant I.V. stimulation have been 
observed in intact animals and man (5-7). Studies show that mild dia- — 
betics man have impairment of only the first or moth ghases of insulin 
secretion (8,9). 

At the time this work was begun there was only one detailea bio- 
logically descriptive "mathematical model" of pancreatic insulin secre- 
(10), although there were several mathematical des- 

‘eriptions in the literature (11-13). That is, for a mathematical des- 


cription to be a biological model it must have been conceptually associated 


with known or predicted biological processes to a reasonable extent. 


Early insulin secretion models were not often biologically very specific 


because they were usually a simplified, small part of an in vivo whole 
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body glucose control model. However, some of the mathematical descriptions 
have been related more specifically to biological processes in later re- 
ports by their authors and thus their differences from biologically des- 
criptive models have become less clear (14,15). | 

we Sen a series of papers, G.M. Grodsky, H.D. Landahl and others 
have developed a model which could account rather well for the multiphasic 
insulin secretion response to an abrupt maintained step-in glucose concen- 
tration (16-20,10). This model was a two-compartment model of the pan- 
creas based on the assumption that insulin was stored in packets which 
varied’ in their secretion threshold to glucose concentration. One com- 
partment* was labile and approximately 2% of a larger stable compartment. 
Other mathematical descriptions of the insulin secretion process include 
the following: 

E. Cerasi produced a lumped variable analog description 
of the in vivo pancreas (11). In this tis there was little attempt to 


correlate the biology to the mathematical functions and the effort could 


be described as rational curve fitting. This single pool model had a rate 


dependent phase. 

3. Porte and Pupo from in vivo studies suggested without a mathe- 
matical formulation a two-pool model for insulin storage and release (5). 

4. Richard O,. Foster presented a iii’ description of the 
in vivo pancreas assuming that there were rate euinitien and concentration 
dependent phases of gucvetion (12). 

Se Srinivasan et al. produced another mathematical description of 
the pancreas with a non-linear rate dependent early phase (13). : 

6. J.R. Guyton in a 1973 thesis produced an elaboration of the 
Foster and Grodsky models (above) (21). 


*For this thesis no distinction is made between the terms compartment 
and pool. 
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7. J. Tiran in his 1974 thesis produced a mathematical version of 
the Porte and Pupo model (above) (22). 

8. E. Cerasi et al. produced a somewhat more biologically descrip- 
tive differential equation (negative feedback) elaboration of the analog 
computer model of Cerasi (above) (14). Other pancreas models or des- 
ortecions exist in the asiepenenie (23-27), but so far the most biologically 
descriptive models are the two pool models of Grodsky (7) and Tiran (22) . 
both mentioned above, 

The work in this thesis was based on the assumption that biologically 
descriptive secretion models are likely to be the most useful in the guid- 
ance of research into the biological mechanisms of secretion. Thus, the 
work that follows will be an attempt to explore biologically descriptive 

alternatives to compartment and pool models both to test the uniqueness 
and quneas applicability of these models and to establish the comparative 
validity of other biological explanations of insulin scanetdens 

In this thesis a two-compartment model has been modified so as to 
better correspond to more recent data developed with the in vitro perfused - 
pancreas, and this modified form has been tested in various experimental 
situations. Experiments which were the most difficult for the two-com- 
partment model to match were used as prime experiments in the construction i: 


of alternate types of models. 


The in vivo distribution and metabolism of insulin in hymans has 
been analyzed, modeled and tested experimentally in order to understand 
the biological significance of the observed secretion kinetics. 

{aeally, the models shown here should not be thought of as a theor- 
etical result, but rather as transient working models, or hypotheses, 
which each summarize our present understanding of a point of view of how 


glucose and stimulated insulin secretion may work. Of course it is quite 
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possible that certain aspects of these models will be retained far into 
the future. However, these models should be thought of primarily as 


an aide to the experimentalist. 
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Il. METHODS 


The in vitro experiments use the perfused pancreas. The technique 
for isolation and perfusion of the rat pancreas has been previously 
described (1,3). The poem upper duodenum, stomach, and spleen | 
were dissected from anesthetized, male, Long-Evans rats. This preparation 
was perfused through the celiac artery with Krebs-Ringer buffer and 4% 
human albumin (furnished by Cutter Laboratories, Berkeley, Calif.) with © 
the aid of a constant temperature and oxygen control system. Perfusion 
rate was 10 ml/min. Glucose, with or without inhibitors, ey ee 
a sidearm syringe. Glucose and drug levels were controlled by mixing 
weighed amounts of these substances into the sidearm syringes and infusing 
them at nis schindnit to give the desired dose. Specimens for immuno- 
reactive insulin assay were obtained from the portal venous effluent. The 
insulin assay is a standard radioimmuno solid phase assay modified from 
the method of Grodsky and Forsham (28). 

6-OH-Dopamine treated rats were given 25 mg/Kg body weight 3 days 


before perfusion to destroy the adrenergic nerve endings in the pancreas 


- as described elsewhere (29). 


The computer work was done in Fortran fon the U.C., SeFe, IBM 360-50 
and the Lawrence Berkeley Laboratory CDC-7600, ana in Focal on the U.C., 
S.F. Department of Biochemistry and Biophysics PDP-12. Most simulation 
plots were done using a printer plotter routine written by Dave Karpinsky 
at the U.C., S.F. computer center. 

The modeling process involved a cyclic interaction between cilia 


and computer simulations. 
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III. Modeling the System 


A. Purpose of Modeling 


All scientific workers are involved in modeling. A hypothesis 
is a model. Thus, any reasonable experiment tests a model which exists 
in the mind of the experimenter. However, if the model becomes 
necessarily complex then the mental image may be less precise than the 
mathematical expression of hs tote. As a working model in a given 
research area becomes complicated, mathematical definition of the model 
together with simulation may be ounieniad for comparison of the model 
to the data. In biomedical research the problem is often complicated : 
by many factors which can be controlled but not eliminated. Thus ,com- 
plex models are frequently required and mathematical modeling becomes 


a useful technique. 
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III. Modeling the System 


B. Techniques of Modeling 


In eth case it was attempted to associate the model as much as 
possible to the biological observations. First a model was conceived. 
Then the model was described mathematically. The following three differ- 
ent approaches were used to determine the values of the model parameters: 
(1) Pure trial and error simulation. In this case the modeler repeatedly 
adjusts model parameters to the best estimated value until a satisfactory 
simulation is achieved. (2) Least square fitting of isolated parts of 


the model to appropriate data. This method usually involves the assumption 


that certain aspects of the data are independent of other aspects. (3) 
Least squares fitting of all model parameters to the data at once. (a) 
Using analytic solutions, (b) using repeated 
solutions, were not possible. This last technique (3b) requres the 
most computer time, however , for non-linear problems, 3a is often not 
possible, 

In each case models are made no more complex than required by the 
data. Thus, if the data is definitive enough a unique solution for the 
values of the model parameters should be achieved. However, most of the 
model solutions presented in this thesis are not unique, and thus equally 


simple alternative models can be constructed and compared. Where the 


models are not unique constraints are listed which are common to all 


models tested. 
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III, Modeling the System 


C. Uniqueness of Biological Models 


When any model is conceived, its uniqueness is a problem in common 
for every worker in every branch of science. In physics and chemistry 
the experiments can often be designed so that a unique model can be 
determined. However, the following factors make non-unique models 
common for the mathematical biologist. (1) The experimental apparatus 
is often provided by nature (in our case the rat pancreas). Thus ,nature 


determines the minimum variability and complexity of the experiments. 


(2) Mathematical techniques are most often used for models which have 
become complex. Complex models are less likely to be unique. 

When a model is not kin it is implied that a group of models 
are possible which could simulate equally Well the available data. To 
eliminate the unsatisfactory models from the non-unique group more de- 
finitive experimental data must be obtained. the similarities and 
these between the non-unique models can guide the experimenter 
in the choice of key experiments. Ideally after enough clever ex- 


perimentation a unique model can be determined. 
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III. Modeling the System 
D. The Models Tested 

General Approach 

After much analysis it was decided that the two most useful opposing 
conceptual views of the pancreatic secretion system were that of a signal 
limited system and that of a storage limited system. For the most part 
a storage limited model existed at the beginning of this work. So that 
comparison of a modified version of the storage limited model to a group 
of signal limited models eventually evolved. 

The following is a brief skematic description of each of the models. 
A detailed mathematical description appears in the —w 

2-Compartment Model 

The 2-compar tment model was by several workers in a series 
of papers, with one of the most recent versions sceeeamnier’ by G. Grodsky 
(7). ‘In that version the model was shown in two inant forms (Fig. 
la & b). Recent on of insulin in- 
dicated that the b form was the more correct (30). The model had been 
tested in the a form, and thus the first step of this study was to write 
and test the b form, where synthesis is directly into the large com- 
partment.” It was found that the b form simulations did not differ 
significantly from the a form simulations. The b form is the version of 
the model which will be called the "storage limited" or "two-compartment 
model" throughout this thesis. 

In the iiciaiiiaiidiasiaas model of the pancreas, insulin is stored in 
a labile compartment, containing 2% of the total insulin. The labile 


compartment is refilled from the larger stable compartment which con- 


? 
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| the small compartment and back respectively. Postulated provisionary 
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tains 98% of the total insulin. Insulin in both compartments is stored 
in packets which vary sigmoidally in their glucose threshold for release. 
The constant m controls the release rate of insulin from the labile 
compartment for packets whose threshold has been exceeded. Paranctere 


Kj(P) and Ky control the rate of movement of insulin from the large to 


parameter P is a function of glucose and time and controls the value of 


K, (P) and in doing so controls the rate of refilling of the labile 
compartment. Thus the model produces the early phase of nibeiatinihe 
secretion by rapid depletion of the small compartment win ei slowly 
rising late phase represents refilling from the large compartment con- 
trolled by the value of provisionary parameter, P. P is assumed to 
reach a citi state value in approximately 1 hour. 

3. The Delta Signal Model | | | 

The delta signal model (Fig. 2a) is an attempt to associate. early 
phase insulin release with the rate of change of a signal as opposed to 
release from a small compartment. The delta signal model is a one com- 
partment model of the pancreas. The compartment is assumed not to be 
significantly depleted during a one hour experiment. The compartment 
volume has been lumped into the constants and thus the insulin secretion 
eatin are represented by the kinetics of the parameter controlling 


release from the compartment. In this case it is assumed that a temporary 


gradient is immediately set up across a structure in the beta cell when 
the cell is presented with a rapid upward step in glucose concentration. 
The temporary gradient is assumed to be the signal for early phase insulin 


release. This model includes a postulated provisionary factor, P, 


| 
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identical to the factor assumed in the two-compartment model as the 
primary controller of late phase secretion. However in this model 
the factor is used to modify the release constant. 

4,  Exciter-Inhibiter (Feedback) Model 

The exciter-inhibiter model (Fig. 2b) is based on the following 
assumptions. A fast upward step of glucose concentration triggers a 
rapid rise of the concentration of an exciter in the cell. The exciter 
(or the glucose) triggers a rise of an inhibiter which rises to a 
concentration approaching that of the exciter. The inhibiter inactivates 
the exciter resulting in the transient early phase release; the rising 
late phase of secretion is controlled by a postulated provisionary fac- 
tor, P, identical to the factor in the beniaiene and the delta signal 
model. 

Mathematically ‘the exciter-inhibiter (feedback) model and the delta 
Signal model differ only slightly. Both are pa compartment Signal . 
limited models which operate thru control of a release parameter. The 
difference between os wnt ae inhibiter is analogous to the size of 
the gradient in the delta signal model. However, the gradient in the 
delta signal model occurs immediately when stimulated by glucose, where- 
as the exciter requires a short but finite time to build up. This mathe- 
matical difference has no significant effect on simulations in this thesis. 
Thus the major difference between the delta signal model and the -exciter- 
inhibiter (feedback) model is conceptual. 

Since the predictions are: equivalent ai} nie limited models simu- 


lations in this thesis are done using the delta signal model, but are 


referred to as "delta-feedback" model simulations. 
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IV. Results Section 


A. Kinetic Glucose Studies < 


‘Introduction’ 

The insulin secretion response of the pancreas to a constant step 
of glucose concentration is multiphasic (1-4) and is therefore a com- 
plex response, In this study it is assumed that by careful examination 
of these and other complex secretion patterns one can gain’ insight into 


the fundamental mechanisms of insulin secretion. The method used is | 


to compare different quantitative models derived from certain kinetic 
control patterns of release and ty them with a 
variety of other stimulatory conditions. 

The basic control pattern consists of constant steps st cneaien 
glucose concentrations continuously perfused for 1 hour (Fig. 3). 
All rate constants, compartment sizes and other parameters are fixed 
from these 7 eet. The test patterns include staircase in- 
fusions (Fig. 4), pulse experiments (Figs. 5,7,&8), ramps (Figs. 12-13) 
and others. The models used will be the two-compartment and the delta— 
feedback models described in the modeling section above. These models 
are thought to be conceptual opposites, i.e. a storage limited model 
versus a signal limited model. It should be stréssed that the storage 
limit mentioned above does not mean a limit in the rate of synthesis 
(which is zoom to be slow) but rather a limit in the movement of insulin | 


stored in the cell to a site available for secretion. 
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Ze Results 
All models were fitted to and produced satisfactory simulations of 
experiments where different glucose concentrations, G, were held con- 


stant for 60 minutes and insulin secretion measured (Fig. 3). That is, 


in the two-compartment model, the size of the small compartment and its 


distribution of threshold sensitive packets were chosen to produce the 


first phase response observed in Figure 3. Parameter values controlling © 


the response of the delta-feedback model to a change in glucose con- 


centration were also chosen to match the experimental first phase behavior 
shown in Figure 3. In both models the second phase behavior was con- 
trolled by the identical iis provisionary factor, P, whose value 
at 60 minutes was matched to the observed secretion rates after 60 minutes 
shown in Figure 3, which were assumed to be steady state values. Thusly, 
P was defined as a slowly responding factor. Without further adjustment 
of any constants in either model, all staircase simulations stint 
experimental error (Fig. 4). For all models the staircase simulations 

. were similar because refilling, which is the major difference between the 


two-compartment model and the delta-feedback model, is not the dominant 


factor in that experiment. In ascending staircase exper inenta, accord- 
ing to the analysis used in the models here, one observes secretion from 
previously unstimulated packets or cells, or because higher glucose levels 
stimulate higher signal levels. 

In the experiments where glucose was infused for 5 minutes at G=100mg% 


and then stopped (G=0) for 5 minutes alternately, only the delta-feedback 


models gave a satisfactory simulation (Fig. 5). In this case the two- 


compartment model with the rate constants developed from data in Fig. 3 


could not be made to refill the labile compartment rapidly enough to account 


| 
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for the large spikes of insulin secretion seen. The delta-feedback model 
simulates this experiment successfully because the capacity to secrete 


a normal first phase of insulin secretion is recovered quickly enough in 


that model. In the delta-feedback model recovery of first phase spiking 
ability is not bound to the value of the postulated provisionary factor, 


P, via refilling as it is in the two-compartment model. That is, in the 


two-compartment model, the rate of refilling of the small compartment 


is controlled by the balance between rate constant from the large to 


the small compartment (a function of P) and the rate of release of 


insulin from the small compartment, whereas in the delta-feedback model 


the rate of recovery of spiking ability is dependent only upon the time 


since the last change in glucose concentration, Effectively the delta- 


feedback model increases the rate of total provision during repeated 


fast increases and rests of glucose concentration beyond the value pre- — 


‘dicted by the two-compartment model. 


To test the negative spike producing capacity of the model an 
experiment was performed and simulations attempted using G=500mg% for 
20 minutes, G=150mg% for 20 minutes, G=500mg% for 10 minutes, and G=150mg% 
for 10 minutes (Fig. 6). The experiment shown in am? was designed 
to produce maximum "refractory" effect or ®negative spike undershoot" in 
_insulin secretion rate obtainable by abruptly reducing the glucose con- 
- centration to a lower though ‘still stimulating concentration. For this 


experiment (and simulations of all other experiments) it was not possible 


to produce a negative spike with the two-compartment model. The reason 


for this is that the release constant of the two-compartment model does 
not change when the glucose concentration is loweredy Initially, the 


insulin secretion rate decreases as the glucose concentration falls below 


/ 
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the glucose threshold for some of the cells (or insulin packets). The 
second effect is the slow movement of the provisionary factor, P, to 
its steady state value for the new glucose concentration. The behavior 


of the release constant may not be changed in the two-compartment model, 


because a constant release rate is necessary to fit the control curves 

where the early phase duration is a constant independent of glucose con- 
centration. The rate of change of P in response to a change in glucose 

level could not be altered in the model because it was fundamentally linked . 

to the slow second phase insulin secretion behavior shown in Figure 3. ’ 


For this experiment, the early phase part of the delta-feedback model — 


elease constant becomes temporarily negative and gives a distinct under— 
shoot at 20 and 50 minutes. This undershoot occurs because the direction 
of the gradient is reversed in the delta signal model, and the inhibiter 
concentration is temporarily greater than the exciter in the exciter-inhibiter 
(feedback) model. 
The experiment also showed a hypersensitization. 
That is the spike seen at 40 minutes in Figure 6 is larger than could be 


accounted for alone by provisionary factor, P. The delta-feedback model 


can not simulate hypersensitization reflected by the spike at 40 minutes, 
because the gradient produced by the glucose change at 40 minutes is 
less, not greater, than at zero time, because the size of the glucose change 


is less and because increased P only provides an increased baseline and 


has no direct effect on spike amplitude. The two-compartment model does 


produce such a large spike due to its property of hyperfilling of the 


labile compartment with high threshold packets during second phase secretion. 


Thus, a modification was made in the delta-feedback model to introduce 


into it second phase hypersensitivity. 
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The modification of the delta-feedback model involved the addition 
of an insulin secretion term which was a function of the product of the 
provisionary factor and the effective gradient caused by a change in. 


glucose concentration. The modified delta-feedback model was then fitted 


to the experiments shown in Figures 3 and 6. 


In further testing both the two-compartment model and the simple delta- 


feedback model produced a satisfactory simulation of the experiment where 


G=300mg% for 60 minutes, G=Omg% for 5 minutes, and then G=300mg% again 


for 5 minutes (Fig. 7). However, the modified delta-feedback model pro- 
duced a spike at 66 minutes which was approximately a factor of three too 


large. This experiment showed that modification of the simple delta- 


feedback model hypersensitized the model's second phase behavior too 


strongly. This result implied an additional complexity in the behavior 

of the pancreas which will be discussed below. | | 
The acdels were then tested against an ahaa where glucose was 

perfused for short periods at G=300mg% and then for short periods at 

G=Omg% alternately (31) (Fig. 8). For this experiment the delta-feedback : 


model gave the most satisfactory simulation. The two-compartment model 


was not able to refill the labile compartment swiftly enough (without a 


modification of the Small compartment size) and the modified delta-feedback 
model showed a too rapidly phase 

With the above results in-nind three new sets of experiments were 
designed. As a control for the effects of glucose concentration in the 
experiment shown in Figure 6 a similar experiment was done with the upper 
G value reduced to G=300mg% (from G=500mg%) . The result of this ex- 


periment was much the same as the experiment shown in Figure 6, i.e. 


only the modified delta-feedback model fit the data in a satisfactory | 
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mannor (Fig. 9). Next the experiment shown in Figure 6, in effect, was 
done replacing the first 20 minute infusion of G=500mg% with G=150mg%. 
The result was a markedly reduced spike at 40 minutes when G was changed 
from 150mg% to 500mg% (Fig. 10). However, after a 10 minute wait at 
G=150mg% another G=500mg% squarewave produced a spike twice as large as 
the previous spike. This led to the hypothesis that "hypersensitivity" 


or “potentiation” of the early phase type of insulin secretion to fast 


changes in glucose concentration may be a process in addition to and/or 


independent of the value of provisionary factor, P. 
The experiment shown in Figure 10 suggests that 10 minutes of ex- 
posure to 500mg% glucose followed by 10 minutes at 150mg% glucose is 


enough to stimulate hypersensitivity, even though the value of P has not 


changed enough to account for the difference in either the two-compartment 


model or the delta-feedback model. Hypersensitivity in the two-compartment 


model depends on the value of the rate parameter controlling the rate of 
flow of insulin from the large to the small compartment which is controlled 
by the.value of P. Hypersensitivity does not really exist in the delta- 


feedback model although the P provides an increased baseline which makes 


spikes of insulin appear larger when P is large. A variable hypersensitivity 
factor, Q(G,t) which responds faster than and is independent of P is a 
possible explanation of increased sensitivity of che pancreas to square- 
waves of 6,9,10, and 11) - a sensitivity 
increase which is dependent upon the pancreas having been exposed recently | 
to a high glucose concentration, That is, the behavior observed in Figures 
6,9,10 and 11 could be explained by a stimulation of a hypersensitivity 


factor which quickly becomes available after exposure and remains available 


for at least twenty minutes. It was notable that the experiment shown 
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in Fig. 10 was matched more poorly by the eines andel simulation 
than were the experiments shown in Figures 8 and 9. In fact in Fig. 10 
the relationship of the two G=500mg% spikes were reversed in the two- 
compartment simulation, i.e. in the simulation the first spike was far a 


larger than the 2nd spike. Thus, the simple provision and potentiation , 


observed in the two-compartment model is not sufficient to account for 
all observed sensitivity increases of the pancreas. For the experiment 
shown in Figure 10 the refilling time is not long enough and the value 
of P not large enough for the model to produce that ties i final spike 
of insulin secretion. For the experiments shown in Figure 10 both the 


delta-feedback model and its modified version did much better than the 


two-compartment model but neither of them gave an entirely satisfactory 


° 


match. 

Finally, a third time control exper inent was done which lent further 
support to the hypothesis that an additional hypersensitivity factor exists 
which may be independent of P and which responds to high glucose concen- 
trations more quickly than P. net time of the initial G=150mg% period 
in annteme 10 was reduced from 40 eo 20 minutes, yet the size of the 


two G=500mg% spikes showed little change in absolute or especially relative 


magnitude (Fig. 11). The simulations matched roughly as before for the 
experiment shown in Figure 10. This supported ne conelusion that the 
~-2nd phase results observed were not so much time dependent as dependent 
upon the type of glucose fluctuations to which the pancreas was subjected. 
Ramp experiments were simulated (Fig. 12 & 13). In the slow and fast 


ramp experiment simulations shown in Figures 12 and 13 respectively it 


could be seen that such ramp experiments did not distinguish well between 


the two-compartment model and the delta-feedback models, because all 
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models gave similar satisfactory simulations for such experiments. In 
the experiment shown in Figure 14 the glucose concentration was increased 
smoothly from G=40mg% to G=100mg% over the first 29 minute period, then 
G was held at 100mg% for 10 more minutes; G was reduced to Omg% for 5 


minutes and finally G was returned to 100mg% for 5 more minutes. The 


two-compartment model simulation showed essentially only low secretion 
at any time during the experiment. Whereas, the delta-feedback model, 
which showed low level activity during the ramp to plateau part of the 
experiment, and then showed a distinct fast spike of insulin secretion 
in response to the final squarewave of glucose. The same type of fast 
spike is observed experimentally for all models when presented without 


a preceeding ramp (Fig. 3). The two-compartment model could not produce 


the type of fast late spike shown in Figure 14 because its labile com- 
partment was depleted by the early ramp to plateau of glucose concentration 


portion of the experiment and this labile compartment could not refill 


rapidly enough at the low provisionary levels that exist within the two- 


compartment model for G=100mg%, It should be emphasized that the ramp | 
7 | 
portion of this experiment did not distinguish between models. Differ- | 
ences seen in simulations during the restimulation were similar to those 


of the other restimulation experiments (Fig. 5 and 8). 


3. Evaluation of the Hypersensitivity Factor 


ta ah-attenct to understand the hypersensitivity factor postulated 
above it was noted that Iverson et al. (32) had ae an increase in 
epinephrine release from the in vitro perfused dog pancreas in the absence 
of glucose. The epinephrine release was inhibited by higher glucose. 


It is known that epinephrine can inhibit insulin secretion and so we 


hypothesize that the epinephrine produced at low glucose concentrations 
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was the reason the hypersensitivity factor did dissipate when the glucose 
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concentration went to zero for 20 minutes. A catecholamine assay was 

not available so the Iverson values for the dog were assumed to be 

correct for the rat. Thus, the ieviet-tiinnin asin in Figures 15,16, and 

17 were performed. Figure 15 is the control experiment somewhat similar 
to the experiment shown in Figure 6 except that it was done with the 
purified perfused pancreas preparation described in section IV.B.2. below. 
Figure 16 shows an experiment which differs from the one shown in Figure 
15 in that glucose concentration drops to 30mg% between 20 and 40 minutes 
instead of Omg% as in the innit, Figure 17 shows an experiment like 

the control experiment except that the rats were pretreated with 6-OH-Dopamine 
to destroy adrenergic nerve endings as described in the anthods section. 
Analysis of the experiments in asin 15,16, and 17 showed no. significant 
differences, suggesting the reported epinephrine release from adrenergic 
nerve endings at low glucose concentrations is not the explanation for 

the dissipation of hypersensitivity phenomena seen above. 


4. Discussion 


Both the two-compartment model and the delta~feedback model (and the 


modified delta-feedback model) matched most of the data, but all of the 


models had difficulty (not always the same difficulty) with certain ex- 
periments, however, each model was helpful in the analysis of the ex- 
perimental results. Table 1 shows a summary of the performance of the 


“Models for each type of experiment. In general it appeared that the 


gradient type of model (delta-feedback model) provided as satisfactory 


a description as the two-compartment model of the behavior of the pan- 


creas. The original delta-feedback model has no potentiation init. A 


simple modification of the delta-feedback model's second phase behavior 


described above is not sufficient to allow complete description of the 
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second phase (potentiation) behavior of the pancreas in all cases. To 
match all of the experimental data it will be necessary to modify the 


second phase behavior of the delta-feedback model in some way so that it 


has the second phase hypersensitivity observed in the in vitro pancreas. 

The above results indicate that, to simulate the hypersensitivity 
correctly, an additional Sector east be included ata on have kinetics 
independent of the slow kinetics of provisionary factor, P. The factor 
must build up within 10 minutes at 500mg% glucose, it must stay up for 
20 minutes, and it must dissipate rapidly when the glucose concentration 
reaches zero or near zero. An alternate explanation might be an increase 
in first phase sensitivity to P under the above conditions. However, 


either change would be an equal complication to the altered model. Char- 


acterization of this behavior is a more complex experimental problem than 


those modeled above and much more experimental data would be required to 


describe it quantitatively in a model. 


The results of simulations by the delta-feedback model including 
negative undershoot in insulin secretion (Fig. 6,9,10, and 11) and in 
simulation of spike secretion of “insulin for low pulses of glucose con- 
centration (Fig. 5 and 14) strongly suggest a sensitivity to rate of 


change of glucose concentration in the mechanism of insulin secretion. 


This will be discussed in detail in a later section. This aspect of 


insulin secretion is consistent with electronic studies of pancreatic 


beta cells (33,34) which show that depolarization of beta cells and the 


number of beta cells initiating action potential discharges increase with © 


increasing glucose concentration. The ion requirements of beta cell 


electrical excitation are similar to-many requirements observed in the 


perfused giant squid axon (35,36). Thus, by analogy if these beta cell 


| 


action potentials are required for first phase spike insulin secretion, 
the disappearance of fast spikes of insulin secretion when the glucose 
concentration is increased slowly may be due to electrical accommodation 
of the peta ‘dels in the population to the depolarization by glucose. 
Conversely ie negative undershoot observed when the glucose concentration 
is dropped quickly to a less effective glucose concentration may be due | 


¢ 


to hyperpolarization of beta cells in the population (37). 


The delta-feedback model is a rate dependence model. It is con- 


sistent with secretion control being accomplished by a rapid flow of 


< 


ions of Catt, Nat, and K* and the: corresponding temporary electrical | 


change. The delta-feedback model is mathematically similar to the two- 


factor model of Rashevsky and Landahl which was found to have many of 
the properties of electrically excitable cells (38,39). However, there 


is little doubt that a sufficiently complex version of the two-compartment 


model could also be matched to the data. (e.g. Guyton in a thesis pro- 


duced negative spikes with a complex 4 compartment alternate to the single 


labile compartment in the two-compartment model used here.) 

Ramp experiments were found to be no test of the two models. This 
is because the threshold distribution hypothesis is essentially cniatunba 
to simple derivative rate dependence, at least for the first stimulation 


before refilling comes into play. # 


- Feedback inhibition is still a possible model for the insulin secretion 
process. The exciter-inhibitor concept shows that this type of model can 
simulate the data as well as the other models. The actual biological © 


| processes represented by the exciter and the inhibitor are yet to be de- 


termined with certainty. Several candidates exist and will be discussed 


in the general discussion section at the end of the results section. 
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All models have the same kind of second phase controlled by the same 
postulated provisionary factor, a mechanism stimulated by glucose but | 
different from that responsible for first phase release. This is consis- 
tent with experimental observations where the first and second phases 
can be differentially stimulated or inhibited by various chemical agents 


(1,40,41,42). 
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IV. Results Section 


B. Rate Dependence 


Ls Introduction 


As noted glucose, when presented to the pancreas as a rapid continuous 


single step, causes a fast but transient release of insulin that terminates 


during the first 5 minutes of the maintained stimulation (Figs. 1-7). 
Subsequently, a second phase occurs, characterized by a slowly increasing 
rate of cenierion, Similar multiphasic responses are also seen in other 
endocrine systems (43,44,45). Some mild diabetics have impairment of 
only the early phase of insulin secretion (8,9), ana Albisser et al. 
have nen it a fast release is necessary for minimizing the insulin | 
required for accurate blood glucose control by a mechanical pancreas (46). 
A sigmoidal dose-response curve is observed when total insulin re- 
leased in either phase is plotted against glucose concentrations used in 
a series of single-step stimulations (10). These types of studies, 
however, do not evaluate the importance of the rate of iii of the 
stimulator as a siena. When glucose is presented as a linearly increas> 
ing ramp, the first phase appears to be smaller but is mre prolonged 
with the decreasing slope of the ramp and ultimately becomes difficult 


to measure in the presence of the gradually increasing second phase 


(Fig. 10). Cerasi (11), Foster (12), and Srinivasan et al. (13) in- 
corporated a rate-of-change signal for first-phase insulin release ina 
pancreatic whole body model, but this was not required in equally effective 


computer models developed from in vitro pancreatic studies (10,17). Al- 


though studies involving comparisons of ramp and step glucose concentra- 


tion changes on first-phase insulin secretion show a clear change of 
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pattern, they have not established whether a fast change in glucose con- 


centration causes more total insulin secretion than a slow change. (This 
is equivalent to asking what has happened to the eurts chan in the slow 
ramp, where it is no longer discernible.) Experimental design is com- 
plicated by the following three special problems: 1) insulin release is 
a non-linear function of glucose concentration; 2) the insulin response 
to glucose produces an overlapping ‘isin Cina and 3) the full response 
of insulin to a given change in glucose ecbeuetvetinn Contin 4to5 
minutes (10). Previous studies relating to the above question could not 
easily control these problems because they were performed in vivo (47) 
or were designed dob cee purposes (48,49). 

2. Methods 

The in vitro perfused rat pancreas preparation with non-recirculating 


phosphate-bicarbonate buffer containing human serum albumin (1%) and 


dextran (4%) were used, as previously described in detail (3). The flow 


rate was 4 ml/min. The experimental design for comparison of glucose 
presentation as a ramp or step configuration is shown in the upper portions 
of panels a and b of Fig. 18. | 

The limits of the change in glucose concentration (from 110 to 180 
mg/dl) were Saabs to restrict studies to the approximately linear region 
of dose-response curve established with this 
preparation (10). | 

Additionally, 110 mg/dl is a good base level because at that glucose 
concentration there is little ieee insulin secretion to complicate 
measurement of first-phase release (10) (see lower portions of 18a and 


18b). 


As previously demonstrated (10,17) and shown here (lower portions of 


18a and 18b), about 5 minutes are required for full insulin secretion at 


v 
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than 5 minutes from the end of each change would show a steadily de- 


__the total time of the change, and the average rate of change are the — 


a maintained glucose step. After that, particularly at high glucose con- 
centrations, late-phase re will interfere. Thus, in comparing both 
cases, insulin secreted for 5 minutes after the end of the change was 

included. | 


It was to be expected that comparing ramp and step periods longer 


creasing percent difference in total secretion and that using less than 


5 minutes would be to the etventen of the ramp because the second step 
of the two-step change would not have had time to give its full early 
phase response. Examination of the data confirmed this. All constraints 
favored either seeing no effect or seeing greater secretion from the | 
ramp, and thus any rate dependence indicated by increased release by | 
the two-step stimulation aay not be the full value. 
Results 
During the experimental period, glucose was increased from the base 
level of 110 mg/dl to 180 mg/dl in two steps (35-mg/dl increments of 5 
minutes duration) or as a ramp (eleven 6.36-mg/dl steps of 1/2 eisute 
duration) . Thus, in compar ing ie tb-step change with the ramp change, 
the limits are in the linear portion of the curve and are the same. Also, 
the average concentration of glucose during the respective experimental 


period period, the total amount of glucose going through the pancreas, 


Cm, : 
same. The only difference is the path (i.e. the size of the step changes 


or the maximum rate of change, which is 5.5 times greater in the two-step 


than in the ramp stimulation) ; the time required to make a concentration 
change was essentially the same for each of the two large glucose steps 


as it was for the eleven small glucose steps. To control the possible 
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change was first, it stimulated secretion of 50% more insulin than the 


why fast changes stimulate more release. 
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influence of sequencing, the two-step stimulation either preceded (panel 
18a) or succeeded (18b) the ramp. A reproducible spike of insulin re- 
lease occurred with each of the two-step stimulations; the ramp experi- 
ments gave no significant spikes but showed a progressive rise of insulin 
of more variable magnitude among experiments. 


As shown in Fig. 19, in the five experiments in which the two-step 


ramp. In the five experiments in which the ramp change was first, the 


two-step change still produced 25% more insulin secretion. When the 


results of all experiments were compared, it was found that 36% more 
insulin was secreted in response to the two-step change than to the eleven- 
step ramp. The data suggest a maximum-rate dependence for glucose-stimulated 
insulin release which, under the conditions of these quper taunts, can 


account for more than 30% of the glucose-related first-phase insulin re- 


lease. 

4. Discussion 

In these experiments it is an ithe to observe directly the 
time of onset of late-phase secretion. However, 5 minutes after the end 
of the change the average secretion rate is actually 23% (not significant) 
higher for the ramp than for the two-step nian. This possible difference 
Suggests that a larger amount of late-phase secretion is not the reason 

We know of no other endocrine system that ne been studied in this 
way. However, rate sensitivity has been discussed in general by several 
suthors (50-53) and specifically with regard to bacterial chemotaxis (54), 
liver adipose cells (55), the pituitary adrenal system (56-58), and the 
pancreas (59). Other biological systems may act similarly: e.g., a 


slowly rising current, much greater than a rapid increase, can depolarize 
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a mew axon without triggering an action potential (60). This type of 

| rate dependence coviires at least two processes which proceed at different 
rates. Among the several types of models that produce this type of be- 
havior qualitatively are the following: 1) a model with a labile com- 
partment that begins to "leak" some of its contents slowly back to a 
seme compartment when secretion is stimulated (and, thus, slow glucose 


concentration changes would cause a greater loss of available insulin 


from the labile compartment than fast changes); 2) a Hodgkin-Huxley type 
of model, which has a fast and slow gate (61); and 3) a model in which a 
rapidly rising exciter and.a more slowly rising inhibitor of secretion 
increase to their respective steady states in response to glucose stimu- 
lation, but with a threshold in the exciter-inhibitor difference required 
for secretion or with secretion dependent on the amount of exciter-inhibiter 
difference raised to a power greater than one. 

The rate-sensor property demonstrated for the pancreas is advantageous 


for effective control systems, and thus may prove characteristic of various 


biologic homeostatic mechanisms. 
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IV. Results Section 


Cc. Diazoxide and DPH 


1. Introduction 

Diphenylhydantoin (DPH) , a drug that reduces the frequency of 
seizures, can also inhibit insulin secretion in man (62-65) and in 
laboratory pancreatic preparations (66-68). Another drug which inhibits 
insulin secretion, diazoxide, has been used purposely to elevate the 
‘blood sugar in certain hypoglycemic en (69-72). At least part of 
this effect of diazoxide is related to inhibition of insulin secretion, 
an effect that has also been demonstrated in vitro (73-75, 68, 76). 

In other studies the responses to DPH ann diazoxide were compared 
during the second phase of glucose-induced insulin secretion (68). It 
was shown that concentratidns of either agent, carefully selected to cause 
a prompt initial inhibition of the second phase of the same magnitude, 
produced characteristically different responses when its perfusion was 
maintained during 20 minutes of a continuous, 60-min, 300 mg/100 ml 
glucose perfusion. DPH caused constant inhibition with return to preper- 
fusion secretion levels immediately after discontinuation of the drug. 
Diazoxide inhibition was characterized by a transient inhibition followed 


by a gradual escape; termination of the drug resulted in a large overshoot 


“Of insulin release. In the compartmental model, assuming a dual action ~~ 


of glucose both on release and on provision of additional insulin (10,17), 
the primary action of DPH was consistent with a cebenittenn inhibition of 
all phases of glucose action, both on release and provision. In contrast, 
the effect of diazoxide appeared to be is ai close to, the final release 


mechanism, with little effect on provision. In a series of experiments we 
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have expanded previous studies to contrast further the effects of these in- 
hibitors, particularly on the initial phase of insulin secretion. In 
addition, we have used the observed stimulus-inhibitor interactions to 
‘Cormulate possible mechanisms of inhibitor effect and to test, with the 
celina the two-compartment hypothesis for insulin secretion. 
2. Results 
Experiment 1. Qualitative effects of inhibitors on early and late 
Secretory phases of glucose-induced insulin secretion are shown in rig. 20. 
Previous studies showed that diazoxide (75yg/ml) and DPH (25yg/ml), infused 
durine the second phase of a glucose stimulation (300 mg/100 ml), produced 
almost complete and similar initial inhibition. Thus, these drug doses 
were used in the present experiments. | 
As shown in Fig. 20, glucose alone caused biphasic insulin vine. 
The simultaneous addition of one of the inhibitors with glucose caused 
their respective ‘drtiein (68) inhibition of the second phase: for 
DPH, second-phase inhibition (about 958) was maintained and secretion re- 
ssetial to normal after termination of the infusion; for diazoxide, escape 
from inhibition occurred within a few minutes and, after termination of 
perfusion, there was a large iain. (i.e., a transient hypersecretion). 
First-phase release persisted, though it was reduced. The following ex- 
_+ ,, periments investigated more closely the effects on the early phase. ; 
Experiment 2. Effects of cf on 
are given in Fig. 21 and Table 2. In a series of separate experiments, 
the effect of preexposure to DPH and diazoxide were compared. The influence 
of each drug on the response to glucose was studied at two different doses 
(25 and 75yg/ml). Drugs were started simultaneously with or preinfused for 


10 minutes before a 5-min infusion of glucose. Results were analyzed by 
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summating insulin values obtained during the first 2-5 min, the summated 
insulin (§ 2-5 immunoreactive insulin (IRI)). 


Control glucose infusions caused typical, early "spike" discharge 


(Fig. 2la). DPH, 25yg/ml, when started taneously with glucose (Fig. 


2lb), inhibited summated insulin. PreinfUSion of the same dose of DPH 
for 10 minutes further inhibited summated early tanelin, compared to in- 
fusions in which DPH and glucose were begun simultaneously (P<.001). 
Raising DPH to 75yg/ml caused no greater inhibition than that achieved by 
DPH preinfused at 25yg/ml. 

Diazoxide, at 25yg/ml and 75ug/ml, inhibited summated early insulin 
secretion (Fig. 2lc). The preinfusion of diazoxide at these doses also in- 
hibited summated first-phast insulin (£ 2-5 IRI). However, the inhibition 
of & 2-5 IRI was no more than with simultaneous infusion of diazoxide, 25yg/ml 
or 75yg/ml respectively. 

— inhibition by low-dose DPH was time and dose dependent, i.e., 
preinfusion produced further inhibition. With diazoxide, however, effects 
were dose dependent, that is, preinfusion of the low dose did not increase 
inhibition. Table 2 statistically summarizes the results of experiment 2. 

| Experiment 3. Influence of inhibitors on response to sequential | 


glucose pulses. In this experiment (Fig. 22 and Table 3) glucose stimulations 


(pulses I and II) were for 7 minutes with two 10-min rest periods, without 


glucose, after each stimulation. Inhibitors were continuously perfused 


10 min before and during the two sequential stimulations. After the second 
rest period, the inhibitor was stopped and glucose alone was given for 10 
minutes (pulse III). 

In control SENN (Fig. 22a), early release was similar in response 


to all three pulses. When either inhibitor was used, there was inhibition 
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of summated release during glucose pulse I compared to control values. 
However, the effects of DPH (Fig. 22b) differed quantitatively from those 

of diazoxide (Fig. 22c). The response to pulse II was less than in pulse I 
with DPH. On the other hand, with diazoxide present, after inhibition of 

the insulin response to glucose pulse I, no further inhibition in summated © 
response occurred during the second pulse. With pulse III (inhibitor removed) 


both inhibitory systems showed recovery, i.e., more insulin was secreted 


compared to pulse II. However, summated early secretion in response to 


glucces gules III was greater in control experiments than when DPH (P<.001) 

or diazoxide (P<.04) were present, indicating that recovery from inhibition 
was incomplete. In addition, the response to pulse III in the DPH-treated 
ein a less than in the diazoxide-treated preparations (P<.01). The 
results of experiment 3 are summarized statistically in Table 3. 

| The studies expand on previous observations that DPH and diazoxide 

have qualitatively of inhibition (68) . These differences 
were evaluated by the two-compartment model previously described and illus- 
trated in Fig. lb. 

As shown in Fig. 23a, the rear hypothesis for the effect of di- 
phenylhydantoin involves reduction by the drug of several components of the 
system. These include terminal secretion (m), the effect of provision (% Be : 
and a reversal of the flow of insulin (increased K') so that the trend is 
for insulin to move toward or remain in the stable compartment.+ 

Fig. 24, a-c, demonstrates the computerized version of the control 
experiments and compares the control simulations with theoretical plots 
of all experiments with inhibitors. The size of the labile compartment 
has been reduced to match the control data early phase. These plots closely 


approximate the experimental observations with DPH: first-phase secretion 
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was partially inhibited (decreased m); preexposure to drug enhanced in- 
hibition (decreased availability of insulin in the labile compartment dus 
to the time-dependent effects of continued decreased x and increased K'); 
and there was no overshoot after discontinuation (though provision was 
stimulated normally, it was ineffective until discontinuation (normal P, 
decreased ¥ )). 

As seen in Fig. 23b, our present hypothesis of the effect of diazoxide 


involves separation of the labile compartment into two compartments. The 


smallest subdivision contains 0.002% of total insulin, and the larger, 1.988%. 


Both have the same qualitative packet threshold distribution. Division of 
the labile compartment was necessitated by an apparent paradoxical difference 
in the effects of diazoxide, depending on whether it was added before or 
during the second phase of secretion. Diazoxide, added during the second 
phase, caused transient inhibition escape and overshoot after termination | 
(indicating a buildup of insulin in the labile compartment by inhibition 
of terminal release (decreased m)) (68). However, in these present 
experiments low drug concentration (25yg/ml) only mildly inhibited the 
release of labile insulin during the first phase of glucose stimulation 

unaffected). Only when diazoxide was increased ‘to 
75ug/ml was reduction of m apparent, inhibiting total early insulin 

(Fig. 21). Response to a third glucose eats after removal of diazoxide 


(Fig. 22, Table 3) showed recovery in sumnated 


tReversal of flow of insulin should not be taken literally, since destruction 
(O<) of insulin from the labile compartment is a mathematically equivalent 
phenomenon and could be a possible effect of DPH. The major suggestion is 
that insulin in the labile compartment is reduced in a fashion which is not 


associated with accumulation of insulin during the inhibitory phase in contrast 


to diazoxide. 


39 
| 
| 


40 


early release. This again suggests involvement of sites closer to the pro- 
visionary phase, with sparing of more distal mechanisms (m relatively un- 
affected). Modification of the theoretical model to include the action of 
diazoxide closely approximated all experimental data (Fig. 24, a-c), assum- 
ing this drug inhibits transfer of insulin from a normally insignificant ; 
intermediate compartment (decreased K,) with other processes remaining 


comparatively unchanged. Thus: first-phase secretion was comparatively 


unchanged (normal to slightly reduced m); second-phase ealeaie was transiently 
inhibited (decreases K,); escape occurred (normal provision to the labile 
compartment with normal x (effect of P) and K); preexposure to drug had 

no effect (no change in availability of insulin in labile compartment, 

normal ¥ and K) ; overshoot appeared after discontinuation of the drug (re- 
lease of previously decreased K,, with normal provision (P), x, and x ). 


Be Discussion 


Based on the two-compartment hypothesis, computer analysis has pre- 
viously approximated a variety of experimental designs related to glucose- 
induced stimuli (10). In the present experiment, we have superimposed the 
inhibitory effects of DPH and diazoxide on the stimulatory effects of glucose. 
Each inhibitor had differential effects on the two phases of insulin secre- 
tion. As seen in experiment 1 and in previous experiments (68), incomplete 
inhibition of the first spike occurred at concentrations of DPH which com- 
pletely inhibited detectable second-phase te secretion. 
even at 75ug/ml, there was distinct early secretion, though this concen- 
tration initially inhibited second phase insulin secretion to almost un- 
detectable levels. 

These data illustrate that the two phases of insulin secretion reflect 


completely different phenomena, a conclusion consistent with that previously 


' 


suggested on the basis of differential stimulation of the two phases by 
glucose and tolbutamide (1) and their differential sensitivity to puromycin 
(3). 

Using data derived from observations of the effects i the inhibitors 
during the late insulin release phase, it was suggested (68) that DPH 
could act by retarding all positive ctamats ete the system, resulting 
in constant inhibition and no overshoot after drug removal. It was also 
suggested that diazoxide inhibited either release from the small compartment 
or a late phase in provision while much of provision continued. Thus, 
initial inhibition behind the blockade produced by the drug would be 
followed by gradual escape as provision increased. After drug removal, an 
overshoot would be observed. Current experiments showing that the first 
phase is only partially reduced by these drugs, at levels by which the 
second phase is almost completely inhibited, provide much more constraint 
on the explanation for the inhibitory effects. 


Our present experiments also suggest that the first and second phases, 


though different in character, may start simultaneously with glucose stimu- 
lation (64). ° Thus, measurements of sumated first-phase insulin release 
made over a finite interval could include a by 
release. This may explain why the height of the first-phase release (mean 
early peak) was less affected by the inhibitors than the summated first-phase 


release. Consideration of this overlapping of the two phases was automatically 


made in the compartmental computer program, ang the computer results closely 


approximated the biological responses. 


The chemical mechanisms of action of diazoxide and DPH are obscure. 


In time it may be possible to relate these chemical actions to the "sites" 


of inhibition indicated by model analysis of kinetic secretion patterns. 
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Diazoxide, among its other actions, affects calcium uptake (77) and cyclic 
AMP levels (78). This drug also causes blockade at specific steps in 
glycolysis (79). DPH is known to affect ion flux, possibly by acting on 
_sodium-potassium-dependent ATPase (80). 

‘Finally, as previously suggested (10), other bases for insulin se- 


cretion models may be applicable, such as feedback inhibition (1) and 


differential signal levels (10). Nevertheless, the compartmental model, 

using a single set of constants derived from a series of glucose dose- 
response curves (10,17), tested with the superimposed effects of inhibitors, 


can simulate a variety of intricate changes in secretion induced during 


experiments of several types of design. 
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IV. Results Section 


D. Insulin Distribution and Metabolism 


Re Introduction 


Until this point, this thesis has involved discussion of models of 


the kinetics and the mechanism of insulin secretion. However, to understand 


the kinetics of the concentration of insulin in the blood in vivo, one 
must combine quantitative tetate for both the secretion of insulin and the 
metabolism and distribution of insulin. Thus, a series of in vivo experi- 
ments were executed and analyzed to evaluate these later phenomena. 
Compartment models for insulin metabolism have been presented by 
several authors (81-83). These models usually consist of 2, 3 or more com- 
partments, one of which corresponds to the blood plasma. A second com- 


| partment may correspond to the extracellular space and a third larger com- 


partment occasionally is included which may correspond roughly to total 


body water. Insulin is thought to distribute among these compartments 

and may decay at different rates in different compartments. There is a dis- 
crepancy between the results in the literature for tracer and non-tracer 
experiments which we believe was resolved by S. Genuth (84). He showed 
that iodine labeled insulin is metabolized at approximately 25% of the rate 
of unlabeled insulin. These metabolic differences combined: with the com- 


‘plex kinetics may have delayed the elucidation of the metabolism and distri- 


bution of insulin. However, one model gave a reasonable description of 
insulin metabolism (83) and our results will be described as they relate 
to that model. 

2. Methods 


A recently diagnosed 27 year old juvenile type diabetic male weight 
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64 Kg. was selected for this study since his endogenous insulin levels were 
relatively fixed at a level below 9 wU/ml and showed no rise with either 
spontaneous or induced hyperglycemia. This facilitates evaluation of the 
distribution and disposal of exogenously administered insulin whose value | 
for half-life measurements and metabolic clearance rate have not yet been 
agreed upon by workers in this area (81,82,8488) . The subject was given 
18 units of insulin on two separate occasions, once by single bolus and once 
by a ateady 4 hour infusion (Fig. 25). Serum insulin was measured at 
intervals by radioimmunoassay using the method of Grodsky & Forsham (28). 
A baseline of endogenous insulin was ented tonne for the patient before 
each experiment. 

3. Results 

Examination of Figure 25 shows that tl/2 of the decay of insulin- 


plasma concentration for the bolus gradually increases with time. However, 


~ 


if one assumes that in the steady infusion in Figure 25, the steady state 
insulin concentration is achieved in three half lives, then the calculated 
t1/2 is smaller chan even during the rapid phase at the start of the bolus. 
Compartment models for insulin sii and metabolism have been presented 
by several authors (81-83). In these experinenta) data we found that linear 
two and three cuastmene eninte with linear destruction rates from all 


compartments like some of those in the literature (81,82), each gave a 


| ~~least square fit with the error per point within the experimental range — 


(<10% error per point) for either the bolus alone or the steady infusion 
alone. 

However, none of the linear compartment models which we tried gave a 
good simultaneous fit for the bolus and the steady infusion data together. 


(All gave >35% error per point). We found that a two compartment model 
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| (Fig. 26a) with a destruction rate of insulin from the blood controlled by 


(Also not shown 7 
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a biomolecular reaction which thus obeys Michaelis-Menten kinetics (Vmax= 
103 yU/min/m1 and Km=418yU/ml) did give a good fit (<10% error per point) 
to all of the data at once (Figure 25). A 2-compartment model with linear 
destruction pita in which an insulin binding protein was assumed to exist 


in the blood instead of a biomolecular reaction controlled destruction 


‘rate did not fit the data better than did the linear 2 and 3 compartment 


models. The model in Figure 26a provides a biologically reasonable ex- 


planation for the wide range of clearance rates and halhf times for insulin 


which now exist in the literature. 
Frost et al. (83) have presented a model simiYar to the one shown in~™ 


Figure 26a. The Frost scom Figure. 26a in that there is both 


a linear and a bimo knsulin from the blood compartment. 


/ 


ics). 
in Fig. 26b to our data 


cretion in 
When we 


(with the same as¢ Bulin secretion used for the 


model in Figure 26a to those of Frost for all para- 
meters except K,, Km a ~und that neither the data of Frost 
nor our data were sufficient to define these parameters uniquely. In fact, 
from the blood compartment to the 2nd compartment still allowed 


moving K, 


a satisfactory fit (Figure 26c). ‘Table 4 shows that these parameters vary 


by a ratio of as little as 4 to l and as much as 12 to l depending on the 


starting points used for the fits. 


Using the steady state equations shown below: 


. £eB (K(B) ) (For our model in Figure 26a) 
I=B_ (K3 + K(B)) (Model in Figure 26b) 
I=B (K, + K(B) - Kp 


a (For our model in Figure, 26c). 


K2 +K3 


ariable endogenous insulin 
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rate did not fit the data better than did the linear 2 and 3 compartment 


(Also not shown, the Frost model allows for variable endogenous insulin se- 


45 


(Fig. 26a) with a destruction rate of insulin from the blood controlled by 
a biomolecular reaction which thus obeys Michaelis-Menten kinetics (Vmax= 
103yU/min/m1 and Km=418yU/ml) did give a good fit (<10% error per point) 
to all of the data at once (Figure 25). A 2-compartment model with linear 
destruction rates in which an insulin binding protein was assumed to exist 


in the blood instead of a biomolecular reaction controlled destruction 


models. The model in Figure 26a provides a biologically reasonable ex- 
planation for the wide range of clearance rates and half times for insulin 


which now exist in the literature. 


Frost et al. (83) have presented a model similar to the one shown in 
Figure 26a. The Frost model differs from Figure 26a in that there is both © 


a linear and a bimolecular destruction of insulin from the blood compartment. 


cretion in normals and, if necessary, in diabetics). 
When we fitted this version of the model in Fig. 26b e our data 

(with the same assumption of a fixed basal insulin secretion used for the 
model in Figure 26a) we got similar values to those of Frost for all para- 
meters except K,, Km and Vmax. We found that neither the data of Frost 
nn our data were sufficient to define these parameters uniquely. In fact, 
moving K, 
a satisfactory fit (Figure 26c). Table 4 shows that these parameters vary 


from the blood compartment to the 2nd compartment still allowed 


by a ratio of as little as 4 to 1 and as much as 12 to l depending on the 
starting points used for the fits. 


Using the steady state equations shown below: 


I=B_ (K(B)) (For our model in Figure 26a) 
I=B (K3 + K(B)) (Model in Figure 26b) 
I=B + K(B) - K, Kp : 
K+K3 ) (For our model in Figure 26c). 
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where I is the steady infusion rate, and B is the steady state insulin con- 
centration of the blood. The metabolic clearance rate (MBC) in ml/Kg/min 


is then equal to I . (Blood Vol.) | 
B (Bd. Wt. in kg.) | 


Each set of values of the parameters obtained by fitting the two models to 
our data (Table 4) were tested against the data from other workers. In 
each case we found the values satisfactory for the data as follows: The 
data of Figures 1 & 3 of Frost et al. (83), the data of Figure 2 of SSnksen 
et al. (87) , meee oo Figures 1,3, and 4 of Sherwin et al. (81). 
4. Discussion 
It is not clear from this study whether the rate limiting bimolecular 
reaction is in the transport, sequestering, or enzymatic breakdown of insulin. 
Isolated liver cell studies in the rat on breakdown of insulin have shown 


a Km about 120 times higher than our estimate of the Km regulating clearance 


concentrations (90). Substrate inhibition (91) and a limitation in trans- 
port or sequestering of insulin have been discussed in the literature (92). 
In one study (83) diabetics were found not to have the saturable bimolecular 


degradation pathway. Further study will be required to determine why our 


patient (n=1) still has a. saturable pathway. It may be because he is a 

new diabetic, more than one saturable ‘cilia exists in normals, or the 

fitted values for model parameters affect this observation. It has been | 
argued (83) the liver of the insulin of 

blood passing thru them, that would not be enough to account for observed 


clearance rates and so there must be other destruction sites. If insulin 


is destroyed at its sites of action, then the model shown in Figure 26c seems 


more conceptually appropriate than the one in Figure 26b or 26a. However, 


| 
(89). However, study of the perfused liver has shown saturation at lower r 
| 


the model shown in Figure 26a is the simplest of the 3 satisfactory models 


shown here and the available data does not seem to be sufficient to determine 


which model is most correct. ‘tt is expected that further esperinente | 
will provide the constraints necessary to distinguish between these 3 models. 
In contrast to another report (88), we find that it is necessary to 
have the results of a steady infusion study of about 30 minutes duration 
in addition to the results of a high concentration bolus study to adequately 
define the insulin binds chicnitertecten of an individual patient. For 
example, one might fit by least squared difference a linear two compartment 
model with only linear destruction rates to data froma high concentration 
bolus experiment alone. In that case, we would expect the resultant model 
to predict too high a steady state blood concentration for an experiment 
where insulin is infused I.V. at a steady rate in the range of e.g. 10 units 


or less per hour. 
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v. General Discussion 


It should be clear to the reader at this point that the mechanism of 
insulin secretion is not well understood. However, the models presented 
here are well defined kinetic hypotheses which must be accepted or rejected 


in the process of elucidating the insulin secretion mechanism. There are 


a large group of processes which might account for the observed secretion 
kinetics. The following is an attempt to list the processes involved in 
association with the appropriate phase of secretion. If the rise of the 
early phase is considered a separate process, then associated processes 
might involve , in some unknown order, the rise of cAMP (93), the influx of 
cat+ (94), the influx of Nat (95), influx of glucose (96), the depolarization 
of the beta cell (95), the initiation of rapid electrical action potentials 


(95), ‘the mobilization of internal bound calcium to a free form (97), the 


inhibition by Nat of mitochondrial ca*t pumping (98), the inhibition of Cat* 
efflux (99), the formation of glucose metabolic products (100), the inter- 


action of glucose and/or a product and/or Catt and/or ATP and/or cAMP 


with any of several receptors (101,100), the suppression of the charge on 

an insulin granule ~ catt (102), the build up of some gradient across some 
structure in the beta cell (electrical or concentration gradient) (41), the 
contraction of microtubules or oe (103), the break down of ATP 

~(40), the binding of vessicles to cell membranes (101), iene of insulin 
oan a labile compartment (10), competition of Catt for sites with polyamines 
(104), movement of granules closer to the membrane (105), increased granule 


velocity (105), decreased cytoplasmic vicosity, and other processes. If 


the fall of the early phase of insulin is a separate process, then it might 


be associated with depletion of a small labile compartment (10) , deactivation 
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of a rapid secretion signal (41), initiation of mitochondrial Ca*t pumping 


(106), binding of cytosolic Catt by the cell tissues (97) , insulin negative 
feedback (107), a change in K*t permeability or other electrical accomo- 
dation (61), a change in Nat permeability (61), local depletion of Catt by 
its loss with insulin (108), receptor desensitization (52), dissipation of 
a gradient across a key cell structure (41), local depletion of cAMP (101), 
atTp (101), glucose (101), or glucose metabolites (101) and other processes. 
The late phase of insulin secretion if a separate process, then it might : 
be due to saturation with or decreased pumping of catt+ by mitochondria (106), 
protein synthesis (3), increased rate of mobilization of stored insulin (10), 
or some other provisionary process. 

No doubt, many of the above processes are involved in secretion but 
at this time it is difficult to eliminate any of these processes from con- 


Sideration. Thus, the function of this thesis has been not to solve, but 


to elucidate the problem. A long series of small atape will be required 
before the above processes can be organized or eliminated with certainty 
or even with high probability. 

It is clear from this thesis that taeah secretion from the pancreas . 
had the following characteristics: 1) the pancreas secretes in a multi- 
phasic pattern in response to a step in glucose concentration, 2) there is 
an early phase with sigmoidal dose dependence, 3) the early phase responsive- 
ness can be increased by previous stimulation (hypersensitization or poten- 
tiation), 4) there is a late phase with sigmoidal dose dependence, 5) the 
hypersensitization of the early phase is a process in addition to and maybe 
independent of provisionary factor P but dependent on many aspects of the 

history of the pancreas, 6) a ramp in glucose concentration presented slowly 


enough gives a monophasic secretion response up to a plateau and 7) a rapid 
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drop in glucose concentration can produce a negative spike in insulin se- 
retion. 
The studies above (see Table 1) show that the two-cotmpartment model 


produces the correct hypersensitivity behavior for only part of the data, the 


negative spike behavior not at all, and has difficulty refilling the small 
compartment fast enough in some cases. However, the delta feedback model 

, has no aypersensstisation behavior and a simple modification was not 
satisfactory. Two points should be clear, 1) either model could be elabor- 
ated to match all of the data included here (but not necessarily — 
experiaente, and 2) the characteristics required of a more appropriate model 
could include properties of both models. 

Thus, a reasonable hypothesis for a better model would be a mixture 

of the two-compartment model and the oitce-taniinaict model (combined model). 


The following is such a model: 


synthesis 50% 50% 


Stable (Q) Labile 


{ (B-1) (P)1* 


This combined model assumes two compartments, but no small compartment, each 


compartment containing 50% of the total pancreatic or cell insulin. The 
early and late phases would be the result of changes in the release con- 
stant of the more labile compartment with the general characteristics of 
the delta-feedback model. Hypersensitivity could aie most of the 
contents of the less labile compartment were shifted into the more labile 


compartment. This process could be controlled by a postulated potentiation 


parameter 0(G,t) via the exchange rate parameter controlling flow of insulin 


from the less to more labile compartment. 
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Due to its larger number of parameters this model would be even tie 
unique than the two-compartment or delta-feedback models. However, given 
enough experimental data it should be possible to test it further. In the 
meantime this model represents a conceptual continuum of models in parameter 
space between the two-compartment and the simple delta-feedback model. For 
as one reduces the less labile compartment in size the combined model be- | 
comes the delta-feedback model. As one reduces the labile compartment 
in size one must hold the release constant more rigid, move the provisionary 
term (P) to act on the refilling constant (K,), to move toward the two- 


compartment model. With enough data one may be able to determine the unique 


values of the parameters in the combined model or show why it is unsatisfac- 


tory. 


The insulin secretion results presented for rate dependence, DPH and dia- 


zoxide show that even the above set of combined delta-feedback and two- 


compartnent models are not sufficient structures in their present form. 

That is, a more appropriate model ake te required to exhibit early phase 

rate dependence, a DPH sensitive refilling, and a third compartment brought ~ 

into being by diazoxide or some equivalent effects. These results represent 

3 additional siiuiiiaiiaaia on tention secretion models and when the process 

becomes well enough understood it will no doubt be possible to produce a | 
model which will be consistent with these constraints also. Once again | 
more detailed data will be required. And once again the above models pro- 

vide a theoretical structure for the experimentalist which can be used to 


form hypotheses and to ask questions about the process around which to plan 


experiments. 


So then how does the above work relate to the controversies of the 


day in the study of secretion: How can these models help the experimentalist? 
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This work provides a holistic view which indicates the minimum complexity 
of the insulin secretion system. It helps to tell experimenters what they're 
looking for (e.g. 2 or more effects of glucose and the list of constraints 
above) and it can help experimenters know where to look. The above model- 
represent a group of structures formed using the set of constraints de- 
termined from the kinetics of insulin secretion. While the kinetics don't 
give a final solution, they do give an outline to the puzzle into which 
other kinds of data can be used to position the pieces. For example, at 
this time the electrical behavior of the beta cell _ been found to exhibit 
multiphasic kinetics in response to glucose which correlate to insulin 
secretion kinetics of the pancreas (109). This means so far as secretion 
and electrical behavior correlate that the techniques for the study of 
excitable cells can be used to study Secretion and that the models _— 


above can be used to understand the electrical response of the beta cell 


to glucose and its relationship to secretion. “There is controversy among 

electrophysiologists as to whether the beta cell action potentials are 

the signal for secretion or the result of secretion by exocytosis. If the 

action potentials are the signal for secretion this might lend support: to 

a signal limited secretion model where the labile compartment is relatively 

large and the kinetics of the release constant actually represent the 

electrical kinetics of the beta cell. | 
However, if the action potentials are the result of secretion,then the 

depletion of a small labile compartment could account for the kinetics of 

the electrical activity. In this way the entese may interact with the 


experiments until secretion is better understood. . 


In addition to trying to understand the mechanism of insulin secretion, 


eventually it is hoped that we may understand the relationship of secretion 
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kinetics to the action of insulin in normals and diabetics. A satisfactory 
insulin clearance model would assist in uriderstanding the insulin secretion- 


action relationship. Thus, we attempted to develop such a model above. 

Modeling in these two lines of research (secretion and clearance) must be 
blended together to fully understand the disease state in vivo. As aeectonea 
above this has. been tried elsewhere with various simplifying assumptions. 

At this time neither secretion nor clearance models have been determined 
which are unique or sufficient. The models in this thesis are based on in 
vitro evidence and may or may not be sufficient to describe other stimulators , 
and inhibitors of insulin release. A whole body insulin secretion model 
would require secretion and clearance modeling of each stimulator and in- 
hibitor of insulin caldese. Since many of these stimulators and inhibitors 
interact with each other, a ‘eingle model including them all would be extreme- 
ly complex. Thus, one would have a large model with a large number of . 
parameters, using mostly data from in vivo studies which are necessarily 
imprecise. The result would be a very non-unique model. This reflects 
science's present true to life experimental view of the system. That is 
the models are limited by vine is known, and thus a complex total body 
model is not appropriate at this time. 

When and if the day aie when effective total body models can be con- 
structed, knowing the ‘iia, a a given individual will require 


excessive testing initially and constant updating as the parameters shift 


with time. Such an effort would probably be worthwhile only if the testing 


were mechanized or for research purposes. The more feasible route for the 
near future is to make appropriately limited models for specific medical . 


problems and to use modeling as a research tool. 
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Because endocrinology, physiology, and in fact all of in vivo biology 


are necessarily complex, properly applied modeling techniques appear to be 


an essential tool for the understanding of these problems and the applica- 


tion of this information to medical treatment. 


4 

4 
| 
3 


Conclusions 


Le Present insulin secretion models must satisfy a complex series 
of events including fast, transient release, second phase, hyper- 
sensitivity and negative response. : 

2. Specific signal limited secretion models work as well as specific 
storage limited models. 

Ze In addition to previously determined constraints (10) an addition- 
al hypersensitivity factor Q(G,t)) (or equivalent elaboration) is 
required for sufficient insulin secretion models. | 
(a) 9(G,t) has faster kinetics than the provisionary factor, P(G,t). 
(b) Release of catecholamines from nerve endings in the in vitro 
perfused pancreas does not account for the fast decay of hyper- 
sensitivity factor (0(G,t) at G=0. 

4. The in vitro seen is sensitive to the rate of danas of glucose 
concentration. 

5. The action of diazoxide, including overshoot behavior when the 
drug is discontinued, can be simulated by a 3* compartment model; 
whereas a two-compartment model is sufficient to describe the action 
of DPH. 

6. Linear compartment models of in vivo insulin distribution and clear- 
ance are not sufficient to simulate both high concentration bolus 
and steady infusion experiments with a single set of values for : 
the constants. 

7. Models with a saturable pathway for the degradation of insulin 
are sufficient to simulate the clearance from the blood of insulin 
presented by high concentration bolus and by steady infusion using 


the single set of values for the constants given in the text above. 


*Storage (i.e. a compar tment) and release of exciter or signal would 
be an analogous way for a signal limited model to produce an over- 
shoot. | 
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No present insulin distribution and clearance model with a saturable 


pathway is uniquely defined by the present data. 


Measurements of insulin levels following both a high concentration 
insulin bolus and a short steady infusion are necessary to adequately — 
define the insulin metabolic characteristics of an individual 

patient. 


A combination of a storage limited model and a signal limited 


insulin secretion model could simulate all of the experiments 


presented in this thesis, but such a model would not be unique. | | 
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Assumptions Common to all Models 


lL. It is valid to use in each model the sigmoidal dose response for the 
amount of early phase glucose stimulated insulin secretion determined from 
the data shown in Figure 3. | 

2. It is valid to use in each model the late phase sigmoid determined from 
the data in Figure 3. 

3. cxswiée is stored in one or more compartment in the beta cell of the | 
pancreas. 

Se A postulated provisionary jon, P(G,t), exists which controls the 

‘late phase secretion rate by changing with first order kinetics to a steady 


value determined from the data shown in Figure 3. 


Notation Common to all Models ——~" ™- 


‘Definitions: 

t = time, t = 0 at time of first non-zero glucose concentration 
| (minutes) . 

G¢t) external glucose concentration (mg%). 

S(G,t) = instantaneous insulin secretion rate (g/min). 

S}(G,t) = one-minute integral insulin secretion rate (jig/min). 

P(G,t) = amount of postulated provisionary factor (ug). 

K = 3.3 constant found by the method of least squares fitting to 
match the sigmoidal relationship between the amount of early 
phase insulin secretion and the glucose concentration ((mg%)). 

“2 = 1.52 x 107 constant found by the method of least squares fitting to 
match the sigmoidal relationship between the amount of wnesy phase 
insulin secretion and the glucose concentration ((mg%)*). 

K = 7.0 constants found ie the method of least square fitting to 


match the sigmoidal relationship between une and glucose 
concentration, G( (mgs) Kp) . 
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2.25 x 101° constants found by the method of least square fitting 


to match the sigmoidal relationship between (P(G,~) and glucose 
concentration, G((mg%) Xp). 


om 0337 the growth and decay rate of the postulated pro- 
visionary factor P(G,t). 


Equations: 
aP(G,t)/dt =~ (P(G,~) -P(G,t)) 
. + GKp) 


Input-Output: 


Input : G(t) 


Output: S(G,t) 


2-Compartment Model 


Assumptions 


Insulin is stored in two compartments. 
Insulin is stored in a small labile compartment from which insulin is 
available for immediate release and in a large stable compartment from 


which the labile compartment is refilled during insulin secretion. 


The rate of refilling of the labile compartment from the stable compartment 


is controlled in part by the value of a postulated provisionary factor, 
P(G,t), which is assumed to follow first order kinetics. 

Insulin is stored in packets in the two compartments. The packets vary 
in their threshold sensitivity to glucose for release. (2(@,0)d@ is 
the amount of insulin in the packets, for threshold level 9® (in units 
of glucose concentration) , between @ and @ + d@. The total amount of 
insulin releasable from the labile compartnant at glucose sini 
G is then: "G 


Z(8,0)de 


Cp 
| 
| 
4. 
X(G,0) = 
@) 
> 
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and X(G,0) is the intial integral distribution function which can be de- 


termined from the data by plotting to total amount of insulin release during 


the early phase against the glucose concentration. 


Synthesis 


98% 2% 
(x(G.1)) 
2-Compartment Model (Contd.) 
Definitions: 
X(G,t) = amount of insulin in the labile compartment as a function of 
glucose concentration G and time t (yg). _ 
Xs (G, t) = amount of insulin in the stable compartment (jg). 
xmax = 1.65 ug, the initial total amount of insulin in the labile 
compartment. 
K, (G,t) = parameter controlling the rate of exchange of insulin from 


the stable to the labile compartment (min -1). 


X(G,0) = initial amount of insulin available in the labile compartment 
as a function of glucose concentration (yg). (i.e. initial * 
integral amount of insulin in packets of threshold G or less). 


Q = glucose threshold of packets of insulin in the labile compartment. 
(mgt) 
z(@,0) = initial distribution of packets of insulin with glucose thresh- 


old for release, @. (yg/mgt) .. 
K,, = .0002 min-!, the rate constant from the stable to the labile 
\ compartment before stimulation with glucose, i.e. when P(G,t)=0.0. 


P(G, ee) 
p 
Precursors 
| 
Vi pf 
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X_ max = 82.5 ug, initial total amount of insulin in stable compartment 


2 (min~+) the effectiveness of the postulated provisionary 
factor, P(G,t). 


1 


A 


.O1 min™*+, the rate constant 


-622 min~1, the release rate constant for insulin secretion 
from the small compartment. 


Equations: 


AK (G,t) /dt = X (G,t) - (Ky + m) X (G,t) 
K, (G,t) K, Pc, t) (G,0) 
X(G,0) = (© 2(0,0)d0 = xmax + G¥) 
0 
(G,0) = X,max + GK) 


dX, (G,t)/dt = Ky X(G,t)-K, (G,t)% (G,t) 


S(G,t) m X(G,t) 


t = 0 
P(G,t) - = 0.0 | | 
(G, t) = (G, 0) 
X(G,t) = X(G,0) 

Delta Signal Model 

Assumptions 

Re Insulin is stored in a single compartment. 
2. The difference in amount of some product of glucose across a key cell 


structure controls the release parameter and thereby determines the rate of 


early phase insulin secretion. Call it structure A. 


| 
| 

j 
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3. The difference across structure A approaches equilibration within 5 
minutes etter a change in extern glucose 
4. There is a Slow changing provisionary component of the insulin se- 


cretion signal which controls the release parameter and which approaches a 
steady state after 60 minutes of constant external glucose concentration. 
oe The transport of insulin to a position in the cell where it is immedi- 


ately available for release is not a rate limiting process. 


G-R = glucose receptor complex 


GI “ oe P(G,t) = provisionary factor 
ucosé 
| Y¥(G,t) = amount of the glucose 
.~product inside of 
structure A 


y'(G) = amount of the glucose 
product outside of 
structure A. 


Nuc. = Nucleus 


Y'(G) is assumed to reach its value and yY' = in 


magnitude. Thus, the magnitude of the difference across structure A is 


Y¥(G,~) -Y(G, t) 


Definitions: 

Y¥(G,t) = amount of a product of glucose on the inside of some cell 
structure (ig). 

Ymax = 1.65 ug, fitting constant chosen before the lease squares 

| fitting of K and C to the sigmoidal relationship between 

¥(G,~) and G. 

, constant determining the rate of increase (or danay) 
of. Y(G,t) inside.of. some cell structure:::. 

D mint, fitting constant relating the difference in Y(G,t) 


across some cell structure to the instantaneous insulin se- 
cretion rate, S (Gt). 
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1.3 (ug/min) ~-, fitting constant relating the product of P(G,t) 
and the limit of Y(G,t) as t goes to infinity, Y(G,~), to the 
instantaneous secretion rate, S (G,t). 


Equations: 

dy(G,t)/dt = D,(¥(G,~)-¥(G,t)) 

Y(G,~ ) = Ymax (note that Y(G,©) is proportional to above.) 
S(G,t) = ((D, 


Initial Conditions: 


t = 0 
Y(G,t) = 0 
P(G,t) = 0 


Modified Delta Signal Model: 


Same at the delta signal model ore except that S$(G,t)=[(D, (¥(G,») -¥(G,t) )+ 
(G,) P(G, t) +D3 (¥(G)-¥(G,t) ) P(G,t) 

DS = 1, fitting constant which adjusts the amount of hyper-— 
sensitivity behavior of the model so as to match the hypersensitization of 


the second spike in the experiment shown in Figure 6. 


Excitation-Inhibition (Feedback) Model 


Assumptions 


1. Insulin is stored in a single compartment. The difference between 
the amounts of an exciter, E(G,t), and Of an inhibitor, I(G,t), molecule at 
some key site inside the beta cell controls the release constant and thereby 
determines the rate of (early phase) insulin were Call it site A. 


2. The amount of exciter at site A quickly approaches a new steady state 


value when the external glucose concentration is changed. 


| 
E 
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3. Several minutes (5) after a fast increase or decrease in the amount of 
exciter the slowly changing amount of the inhibitor approaches the amount 
of the exciter.. 

4. There is a very slowly changing, internally mobilized, resultant of 
the exciter to which the inhibitor concentration does not respond which 
controls the release parameters and which approaches a steady state after 

60 minutes of constant external glucose concentration. 

5. The transport of insulin to a position in the cell where it is immedi- 


ately available for release is not a rate limiting process. 


‘G-R = glucose receptor complex 
E = exciter 
I = inhibitor 
P = provisionary factor 
Definitions: 
E(G,t) = amount of exciter at nei sienna G and time e (ug) . 
1(G,t) _ = amount of inhibitor at glucose concentration G and time t (ug). 
Emax = 1.65 ug, constant associating the maximum integral amount of 
insulin in the early phase to the upper limit of the sia ceed 
relating E(G,~).and G. .; 
A, | = 3.0 min- 1. constant determing the rate of increase (or decay) 


of the amount of exciter, E(G,t). oi 


Gluc 
| 


H = .6 min-!, constant determining the rate of increase (or decay) 
1 of the amount of inhibitor, I(G,t). 


H = .24 (g/min) ~!, fitting constant relating the product of 
postulated provisionary factor P(G,t) and E(G,t) to the in- 
stantaneous secretion rate S(G,t). 


H = 1.0 min -1 fitting constant relating the difference between . 
the exciter E(G,t) and the inhibitor I(G,t) to the instantaneous 
insulin secretion rate S(G,t). 


N = .95, constant determining the minimum secretion rate at the 
end of the early phase of secretion (no units). 


Equations: 
dE(G,t)/dt = H, (E (G,~)-E(G, t) ) 


(G,t) /at = 

E(G,”) = Emax (c+6* ), (note that is proportional to X(G,0) 
of the 2-compartment model. ) 

I(G,>) =N E(G,~) 

s(G,t) (E(G,t)-1(G,t)) +H,P(G, t)E(G,t)) 


Initial Conditions: 


t = 0 


| 

| 

E(G,t) = 0 

I(G,t) = 0 | 

P(G,t) =0 

| 

| 

| 

| 

| 
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Legends 
FIG. l. Theoretical scheme for insulin secretion incorporating the 


concept of a distribution of insulin units or packets which 


release insulin when their threshold to glucose is reached. 


Open circles in the small compartment represent the distri- 
bution of the threshold-sensitive packets. Broken tists | 
in the large ‘compartment are for artistic balance; the 
distribution characteristics of this compartment are un- 
known. Figures la and b show the form of the model if new 


insulin preferentially enters the small (figure a) or the 


large (figure b) compartments. (From Grodsky (10)). 


FIG. 2(a). Theoretical scheme for insulin secretion incorporating the 
concept”that insulin is stored in a single compartment and 

_ that the release rate is controlled by the kinetics of a 
Signal. The early phase behavior of the signal is controlled  -. 
by the difference in the amount of some product of glucose 
across a key cell structure. The late phase part of ne 
‘Signal is controlled by the value of the slowly changing 
provisionary factor P. (b) A scheme for insulin secretion 
Similar to that in 2(a) (above) except that the early phase 
part of the release parameter behavior is controlled by the 
difference in the amount of exciter and inhibiter at the 


key sites. The exciter is produced rapidly in presence of 


glucose while the inhibiter is either stimulated directly 
by glucose or by feedback through the exciter. 


| 

| 

| 


UMI 


FIG. 3. 


Insulin secretion during single-step continuous perfusion 
with glucose. Solid lines show experimental results. Broken 


line shows results of computer simulations with the 2-com- 


partment model. Dashed line shows simulation with the delta- 


FIG. 4(a). 


FIG. 5(a). 


feedback model (the modified delta signal model did not vary 
from the dashed line). ‘The line of alternating dots and 
dashes shows the results of simulation with the exciter- 


inhibiter model. 


The experimental insulin‘ secretion rate plotted against time. 
Insulin secretion during staircase stimulations with 50 mg/dl 
increments of glucose. § is mean + S.E., n=7. (b) Two- 
compartment model computer simulation of exper iment shown in A. 
(c) Delta-feedeck model simulation of the experiment shown in 
A. (da) Modified delta-feedback ota simulation of the experi- 
shown in A. | 
The experimental secretion rate plotted against time. Effect 
of 4 sequential 100 mg/dl glucose Square waves of 5 minutes 
uration with 5 minute periods of 0 mg/dl glucose between. §& 
is mean + S.E., n=3. (b) Two-compartment model computer simu- 
lation of the experiment shown in A. (c) Delta-feedback 


model simulation of the experiment shown in A. (d) Modi- 


_ fied delta-feedback model simulation of the experiment shown 


FIG. 6(a). 


in A. 


The experimental insulin secretion rate plotted against 


time. Effect of two sequential 500 mg/dl glucBée stimulations. 
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_ FIG. 7(a). 


FIG. 8(a). 
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Comparing the summated ooveetion in minutes 4 through 8 of 

the first peak to the summated secretion in minutes 4] ehrough 
45 of the second peak gives average vidiaiie of 84.6 ug and | 
132.8 ug beaeiectivete. P is less than 02 using paired sta- 
computer simulation of the experiment shown in A. (c) Delta- 
feedback model simulation of the exper iment shown in A. 

(a) Modified delta-feedback model simulation of the experi- 


ment shown in A. 


The experimental insulin secretion rate plotted against time. 
Insulin secretion during a’ 300 mg /al glucose step followed by . 
a 5 minute rest at 0 mg/dl followed by 5 minutes at 300, mg/dl. 

£ is mean + S.E., n=9. (b) Two-compartment model computer simu- 
Oe of the experiment shown in A. (c) Delta-feedback model 
simulation of the experiment shown in AL (da) Modified delta- 


feedback model simulation of the experiment shown in A. 


The experimental insulin secretion rate plotted against time. ; 
Insulin secretion during two 7 minute 300 mg/dl glucose square 
waves followed by 10 minutes rest at 0 mg/dl followed by an 

minute 300 mg/dl glucose square wave. the mean + S.E. 


=6. (b) Two-compartment model computer simulation of the ex- — 


‘periment shown in A. (c) Delta-feedback model simulation of the 


experiment shown in A. (da). Modified delta-feedback model simu- 


lation of the experiment shown in A. 


‘ 

| 
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FIG. 9(a). 


FIG.10(a). 


>. 
The experimental insulin secretion rate plotted against time.- 


Effect of two sequential 300 mg/dl glucose stimulations. 
Comparing the summated secretion in minutes 4 through 8 of 

the first peak to the summated secretion  ainwie 41 through 
45 of the second peak gives average values of 66.5 ug and 
106.6 ug respectively. P is less than .1l paired 
atatietics. 6.8., n=4. (b) Two-compartment model 
computer simulation of the experiment shown in A. (c) Delta-— 


feedback model simulation of the experiment shown in A. (qd) 


Modified delta-feedback model simulation of the exper iment 


shown in A. 


The experimental insulin secretion rate plotted against time. 
Effect of two 10 minute sequential 500'mg/dl glucose square 
waves following 40 minutes at 150 mg/dl glucose and separated 


by a 10 minute rest at 150 mg/dl glucose. Comparing the summa- 


. ted secretion during the peaks in minutes 41 through 45 to 


FIG.11(a). 


that in minutes 61 through 65 gives average values of 63.0 ug 
and 101.6 ug respectively. P is less than .0l1 using paired 
Statistics. is mean + S.E., (b) Two-compartment model 
computer simulation of the experiment shown in A. (c) Delta- 
feedback model simulation of the experiment shown in A. (d) 
Modified delta-feedback model a of the experiment 


shown in A. 


The experimental insulin secretion rate plotted against time. — 
Effect of two 10 minute sequential 500 mg/dl glucose square 


waves following 20 minutes at 150 mg/dl glucose and separated | 
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‘by a 10 minute rest at 150 mg/dl glucose. Comparing the sum-_ 


FIG.12(a). 


FIG.13 (a). 


-FIG.14(a). 


-Delta-feedback model simulation of the experiment shown in A. 


mated secretion during the peaks in minutes 21 through 25 

to that in minutes 41 through 45 gives average values of 82.0 

and 159.3 respectively with P less than .005 using paired statistics. £ 
is mean + S.E., n=5. (b) Two-compartment model computer simu- | 

lation of the experiment shown in A. (c) Delta-feedback model 
simulation of the experiment. shown in A. (d) Modified delta- 


feedback model simulation of the experiment shown in A. 


The experimental insulin secretion rate plotted against time. 
Effect of glucose administration as a slow ramp. Slope of 
about 5 mg/dl/min. & is mean + S.E., n=4.  (b) . Two-compartment 
model computer simulation of the experiment shown in A. (c) 
Delta-feedfack model simulation of the experiment shown in A. 
(a) Modified delta-feedback model simulation of the experiment 
shown in A. | : | 


- 


The experimental insulin secretion rate plotted against time. 
Effect of glucose administration as a rapid ramp. Slope of 
about 50 mg/dl/min. £ is mean + S.E., n=3. (b) Two-compartment 


model computer simulation of the experiment shown in A. (c) 


(d) Modified delta-feedback model simulation of the experiment 


shown in A. 


The experimental insulin secretion rate plotted against time. 
Effect of glucose administered as a slow ramp from 40 to 100 
mg/dl followed by a 5 minute rest followed by a square wave of 


100 mg/dl glucose. & is mean + S.E., n=7. (b) Two-compartment 
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model computer simulation of the experiment shown in A. (c) , 
Delta-feedback model simulation of the experiment shown in A. 


(d) Modified delta-feedback model simulation of the experiment 


shown in A. 


FIG.15 (a). The experimental imsulin secretion rate plotted against time. 


Effect of 20 minutes at 500 mg/dl glucose followed by 20 minutes 


at 0 mg/dl glucose followed by two 10 minute:square waves of 


500 mg/dl glucose separated by a 10 minute rest at 0 mg/dl 
glucose. Summations of insulin secretion in early phase and 
did not differ significantly. 3% is mean + S.E., n=5. (b)_ 
‘Two-compartment model computer atmmlation of the experiment 
shown in A. (c) Delta-feedback model simulation of the ex- 
periment shown in A. (d)} Modified delta-feedback model simu- 


lation of the experiment shown in A. 


FIG.16(a). The experimental insulin secretion rate plotted against time. 


Effect, of 20 minutes at 500 mg/dl glucose followed by 20 min- | | 


utes at 30 mg/dl glucose followed by two 10 minute square 
waves of 500 mg/dl glucose separated by a 10 minute rest of 
30 mg/dl epee. é is mean + S.E., n=5. (b) Two-compartment 
model computes simulation of the experiment shown in A. (c) | 
‘Delta-feedback model attaches of the experiment shown in A. 
(a) Modified €elta=feedback model simulation of the experiment | 
shown in A. | 
FIG.17(a). The experimental insulin secretion rate plotted against time. 


Effect of 20 minutes at 500 mg/dl glucose followed by 20 min- 


| 
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FIG.18(a). 


FIG. 19. 


eten at 0 mg/dl giuscee followed by two 10 minute square waves 
of 500 mg/dl glucose separated by a 10 minute rest at 0 mg/dl 
glucose. Summations of insulin séoretion in early phase and 
did not differ significantly. Eis mean + S.E., n=4. The 
rats were treated with 6-OH-Dopamine 3 days before perfusion 
of the pancreas. (b) Two-compartment model computer simulation ~ 
of the experiment shown in A. (c) Delta-feedback model simu- 
lation of the experiment shown in A. (d) Modified delta- 


feedback model simulation of the experiment shown in A. 


The pattern of glucose perfusion and the plot of insulin 
secreted by the pancreas as functions of time. A two-step 
simulation of identical limits and duration was presented 
first and an eleven-step, ramp-type stimulation was presented 
second. (b) A plot of the alternate experimental design in 


which the order of the stimulations was reversed. 


Total insulin secreted when glucose was presented at different | 


maximum rates of change. Data for each bar were derived from 


the insulin secreted during the shaded periods shown left to 
right in Fig. 18(a) and (b) and the combined amounts of Fig.18 
(a) and (b). Note that in each case more insulin is secreted 
in response to a two-step change than in response to an eleven- 
step, ramp-type change. Probability values were calculated 


using paired Student's t tests. 
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Table 1. Comparison of the performance of each model in simulation 
of the various categories of experimental kinetic insulin 

secretion patterns. + indicates a satisfactory simulation. 


- indicates an unsatisfactory simulation. 


Fig. 20. Effects of inhibitiors on biphasic glucose induced insulin 
release. In this and subsequent figures DPH=diphenylhydantoin. 


Shaded bars represent glucose (300 mg/dl)... 
| Fig. 21. Effect of inhibitors on early release of instlin. 


Fig. 22. Influence of inhibitors on insulin response to sequential 


glucose pulses. Shaded bars represent glucose (300 mg/dl) . 


Fig. 23. | Two-compartment model for glucose-induced insulin secretion. 
Percentages beside symbols in DPH and diazoxide (8) models 
indicate postulated percentage alteration of control values 


(fig. lb) ‘by inhibitors. 


Fig. 24. (a) Computerized simulations of biological data from fig. 20, 
based on 2-compartment model . (b) 
of biological data from fig. 21, based on 2-compartment 
model. (c) computerized simulations of biological data from | “aie | 


fig. 22, based on 2-compartment model. 


Table 2. Effect of prior infusion and dose of drugs on inhibition 
of early insulin release. Values are means + S.E. Insulin, 


ng/ml. DPH=diphenylhydantoin. SIM=inhibitor and glucose started 


Simultaneously. PRE=inhibitor infused for 10 min prior to 


| 
| 


Table 
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coinfusion with glucose. In all experiments, glucose 300 


mg/dl was the stimulus. £2-5=mean sum of values in 


the first 2-5 min after the start of glucose infusion. 


Insulin response to sequential pulses of glucose (300 
mg/dl). Values are insulin, ng/ml + S.E. NS=not significant 
(P>.05). . Summated=mean sum of values in the first 2-5 min 


after the start of glucose pulse. See text and Fig. 22 for 


description of experimental design. 
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Insulin blood concentrations from a single recently diagnosed 
diabetic patient. Q@=values for an 18 U I.V. bolus. given at | 
t=0, and the solid line is the computer fit using the model 


shown in Figure 26. A=values for a ateany 4 Zs V.. infusion of 


insulin at a rate of OF U/Kg/4 hr., and the dashed line is 


the computer fit using the model shown in Figure 25. The 


patient's weight was 64 Kg. 


A model of insulin distribution and degradation in humans. 


For the steady I.V. insulin infusion experiment the infusion 
rate, I=.07 U/Kg/4 hr. The patient weighed 64 Kg. For the 
I.V. bolus experiment I=0, and the initial blood insulin 
concentration, B(t=0) was assumed to be the 18 U dose divided 
by the blood volume. 

(a) For our patient K,=.044 min., K5=.027 min., vmax=103 


ywU/min/ml, and Km=418 yU/ml. If passive transport is assumed 


_ then the least squares fit gives a plasma associated volume 


of 7.07 liters, a 2nd compartment of 11.39 liters, and an 
exchange coefficient of .62. The compartment 2 volume is not 
incompatible with that of the extra-cellular space. 


(b) An insulin metabolism model similar to a model by Frost 


et al. (9). The values of the parameters depended on the 


starting values of the fit. See Table 4. 
(c) An insulin metabolism model with linear iets 
at a site other than the liver or kidney. The values of 


the parameters depended on the starting values of the fit. 


- See Table 4. 


q 
{ 


Table 4. Comparison of parameter values for insulin distribution 
| and clearance models fitted to data from one previously 


undiagnosed diabetic. 
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Table 


Conditions 


Control 


25 


75 


 Diazoxide 25 


75 


Control 
“SIM 
PRE 


II 


<.001 
<.001 


<.001 
<.001 


[+ 1+ 
N 


<.O1 
<.02 
<.001 
<.001 


=) 


P (Paired Analysis) 


II III III 


Control 
(n=6) 


DPH,25 yg/ml 


(n=5) 


Diazoxide,75 


ug/ml 
(n=6) 


105+18 


3+1 


111416 NS NS NS 


2147. .03 NS <.03 


65412 NS NS 


109 
| 
Drug Dose a n Summated P vs. 
ug/ml (= 2-5) 
20 83 +4 
DPH 5 33 + 4 
2 
SIM 3 19 
PRE 3 13 | 
pS SIM 6 57 
PRE 6 52 | 
| SIM 6 39 
PRE 12 38 
| 
Table 3. 
| | 
| 
Summatea 90+3 | 
Summated 1243 
Summated 41+4 36+4 


Bolus Data | 
A Steady Infusion Data 
— Bolus Computer Fit 
--- Steady Infusion Computer Fit 
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FiG. 25 ° 
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| 
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| lll 
Vmax 
Km +B 

KI | 
K(B) 

| 

I 

K3-.1K(B) 

KI 

KB) K3 


‘TABLE 4 
| = COMPARISON OF MODEL PARAMETER VALUES 
4 | Liters 
Min-1 Min-1 Min-1 pu/ml1 wU/min/ml Blood Re- LSD Vmax 
K> Km Vmax lated Vol.. ‘ Km 
| Model in Fig. 26a  .0269 0.0 417.78 103.45 7.071 2245 
Frost et al. orig- ve 
normals) : .058 -109 44.3 16.15 N.A. 365 
Frost et al.orig- 
inal values (for | 
| diabetics) “0506  #£.0619 --- --- 9.363 N.A. 
Model in Fig. 26b (using 
Frost starting values, 
for normals) 052 0248 61.1 22.9 
in Fig. 26b 
(using 26a start 
values) 043 371 | 88°. 7 6.82 0990 239 
Model in 26b (using | | | 
intermediate start- 3 | 
ing point) 0456 .0246 -0373 183.32 46.8 6.33 
‘Model in 26c (using 26a | 
«066 -021 -008 339 87.7 | 6.24 .0953 6259 
Model in 26c ( using 
Frost starting 
normals) -0892 <0122 . 150.6 38.6 6.33 -0975 256 
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